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ABSTRACT 

This  study  analyzes  the  flight-level  data  collected  by  research  aircraft  that  penetrated  the  eyewalls  of 
category  5  Hurricane  Hugo  (1989)  and  category  4  Hurricane  Allen  (1980)  between  1  km  and  the  sea  surface. 
Estimates  of  turbulent  momentum  flux,  turbulent  kinetic  energy  (TKE),  and  vertical  eddy  diffusivity  are 
obtained  before  and  during  the  eyewall  penetrations.  Spatial  scales  of  turbulent  eddies  are  determined 
through  a  spectral  analysis.  The  turbulence  parameters  estimated  for  the  eyewall  penetration  leg  are  found  to 
be  nearly  an  order  of  magnitude  larger  than  those  for  the  leg  outside  the  eyewall  at  similar  altitudes.  In  the 
low-level  intense  eyewall  region,  the  horizontal  length  scale  of  the  dominant  turbulent  eddies  is  found  to  be 
between  500  and  3000  m,  and  the  corresponding  vertical  length  scale  is  approximately  100  m.  The  results 
suggest  also  that  it  is  unwise  to  include  eyewall  vorticity  maxima  (EVM)  in  the  turbulence  parameter  esti¬ 
mation  because  the  EVMs  are  likely  to  be  quasi-two-dimensional  vortex  structures  that  are  embedded  within 
the  three-dimensional  turbulence  on  the  inside  edge  of  the  eyewall. 

This  study  is  a  first  attempt  at  estimating  the  characteristics  of  turbulent  flow  in  the  low-level  troposphere  of 
an  intense  eyewall  using  in  situ  aircraft  observations.  The  authors  believe  that  the  results  can  offer  useful 
guidance  in  numerical  weather  prediction  efforts  aimed  at  improving  the  forecast  of  hurricane  intensity. 

Because  of  the  small  sample  size  analyzed  in  this  study,  further  analyses  of  the  turbulent  characteristics  in  the 
high-wind  region  of  hurricanes  are  imperative. 


1.  Introduction 

Turbulent  transport  processes  are  believed  to  play  an 
important  role  in  the  intensification  and  maintenance 
of  a  hurricane  vortex  (e.g.,  Emanuel  1995, 1997;  Persing 
and  Montgomery  2003;  Smith  et  al.  2008;  Bryan  and 
Rotunno  2009a;  Rotunno  et  al.  2009).  However,  because 
of  safety  constraints,  direct  measurements  of  turbulence 
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in  the  inner-core  region  of  intense  hurricanes  have  been 
essentially  nonexistent  and  observational  efforts  have 
been  focused  almost  exclusively  outside  the  eyewall  re¬ 
gion.  Turbulence  measurements  in  the  high-wind  region 
of  the  storm  would  enable  a  proper  assessment  of  subgrid- 
scale  parameterizations,  which  include  boundary  layer 
parameterizations,  used  in  numerical  models  for  hurri¬ 
cane  intensity  prediction  (e.g.,  Braun  and  Tao  2000; 
Nolan  et  al.  2009a, b;  Smith  and  Thomsen  2010). 

Moss  and  Merceret  (1976,  1977)  and  Moss  (1978)  de¬ 
scribed  one  stepped-descent  flight  pattern  in  the  bound¬ 
ary  layer  of  the  periphery  of  Tropical  Storm  Eloise.  This 
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research  flight  measured  momentum  fluxes  at  different 
levels.  Almost  30  yr  later,  the  Coupled  Boundary  Layer 
Air-Sea  Transfer  (CBLAST)  Experiment  conducted  dur¬ 
ing  the  2002-04  hurricane  seasons  provided  valuable 
data  that  contain  measurements  of  turbulent  fluxes  of 
both  momentum  and  enthalpy  in  the  hurricane  bound¬ 
ary  layer  between  the  outer  rainbands  (e.g.,  Black  et  al. 
2007;  Drennan  et  al.  2007;  French  et  al.  2007;  Zhang  et  al. 
2008).  However,  all  of  these  boundary  layer  turbulence 
observations  have  been  restricted  to  surface  wind  speeds 
<30  ms-1  and  to  regions  well  outside  the  eye  wall. 

Using  the  data  collected  by  a  National  Oceanic  and 
Atmospheric  Administration  (NO  A  A)  WP-3D  research 
aircraft  (N42RF  hereafter)  during  the  eyewall  penetra¬ 
tion  of  category  5  Hurricane  Hugo  (1989)  at  an  altitude 
of  around  450  m,  Marks  et  al.  (2008)  recently  presented 
the  first  in  situ  observations  of  the  structure  and  evolu¬ 
tion  of  the  coherent  vortex  features  in  the  eyewall,  such 
as  the  eyewall  vorticity  maximum  (EVM),  which  was 
observed  to  be  on  the  inner  edge  of  the  eyewall  re¬ 
flectivity  maximum.  During  the  eyewall  penetration, 
N42RF  was  flying  within  the  inflow  layer  according  to 
the  Doppler  radar  observations. 

The  flight-level  data  collected  during  this  Hugo  mis¬ 
sion  are  among  the  very  few  measured  in  the  low-level 
troposphere  of  an  intense  hurricane  near  and  inside  an 
intense  eyewall.  Another  observation  of  this  kind  was  the 
flight-level  data  collected  on  6  August  1980  in  Hurricane 
Allen  (1980).  The  data  from  these  two  flights  provide 
a  unique  opportunity  to  estimate  the  characteristics  of 
the  turbulence  near  and  underneath  the  eyewall  of  an 
intense  storm.  One  objective  of  this  study  is  to  esti¬ 
mate  turbulence  parameters  such  as  momentum  flux, 
turbulent  kinetic  energy  (TKE),  and  eddy  diffusivity  in 
the  high-wind  regime  whose  surface  wind  speeds  exceed 
30  m  s”1.  Another  objective  is  to  estimate  the  spatial 
scales  of  the  energy-containing  turbulent  eddies  in  the 
low-level  eyewall  using  detailed  spectral  analysis. 

An  outline  of  the  remaining  sections  of  this  paper  is  as 
follows.  In  section  2,  we  give  a  brief  description  of  data 
and  the  analysis  methodology.  Section  3  describes  the 
potential  errors  involved  in  the  data  analysis  and  the 
method  employed  herein  to  assure  quality  control.  In 
section  4,  we  present  the  results  of  the  spectral  analysis. 
Section  5  presents  the  results  of  the  estimated  turbu¬ 
lence  parameters  and  a  detailed  error  analyses.  This  is 
followed  by  section  6,  which  summarizes  the  main  find¬ 
ings  and  the  limitations  of  the  results. 

2.  Data  and  analysis  method 

As  mentioned  earlier,  the  data  used  in  this  study  are 
from  two  research  flights,  one  into  category  4  Hurricane 


Allen  (1980)  and  the  other  in  category  5  Hurricane  Hugo 
(1989).  We  analyzed  the  flight-level  data  from  the  period 
of  missions  before  and  during  the  eyewall  penetrations 
when  N42RF  was  flown  at  nearly  constant  radar  altitudes 
below  1  km.  Wind  velocity  data  were  corrected  for  air¬ 
craft  motion,  measured  with  an  Inertial  Navigation  System 
(INS),  and  a  Global  Positioning  System  (GPS).  Surface 
(10  m)  wind  speeds  are  estimated  using  measurements 
from  a  nadir-pointing  stepped  frequency  microwave  ra¬ 
diometer  (SFMR;  Uhlhorn  et  al.  2007).  Note  that  there 
is  no  SFMR  data  during  the  flight  mission  into  Hurri¬ 
cane  Allen  because  SFMR  was  not  installed  at  that  time. 

The  data  from  the  flight  in  Hurricane  Hugo  on 
15  August  1989  have  been  published  recently  by  Marks 
et  al.  (2008).  That  study  provided  a  detailed  description 
of  the  synoptic  conditions,  storm  structure,  the  experi¬ 
ment,  and  data  quality,  etc.  In  this  work,  we  used  the 
data  from  the  portion  of  the  mission  between  1720  and 
1729  UTC.  Flight-level  parameters  of  interest  include 
aircraft  altitude  and  wind  speed  in  the  three  Cardinal 
directions  (Fig.  1).  Also  shown  in  Fig.  1  is  the  surface 
wind  speed  estimated  by  the  SFMR.  The  gray  lines  in 
each  panel  of  Fig.  1  represent  the  three  time  intervals 
investigated  here.  The  first  time  interval  was  taken  out¬ 
side  the  eyewall  and  was  chosen  to  compare  against  the 
results  from  CBLAST.  The  second  time  interval  corre¬ 
sponds  to  the  eyewall  penetration  up  to  the  inner  edge  of 
the  eyewall.  Although  this  flight  segment  includes  the 
rapid  increase  in  the  tangential  wind  with  decreasing 
radius,  it  excludes  the  EVM.  The  third  time  interval 
covers  both  the  eyewall  penetration  and  the  EVM.  Time 
intervals  1  and  2  are  judged  suitable  for  turbulent  flux 
calculations  according  to  the  detailed  spectral  analysis  that 
will  be  discussed  in  sections  3  and  4.  Time  interval  3  on 
the  other  hand  is  not  suitable  for  flux  calculation,  but  it  is 
nonetheless  used  for  comparison  purposes  to  estimate 
how  the  EVM  may  affect  the  three-dimensional  turbu¬ 
lence  characteristics. 

The  data  used  in  this  study  from  the  flight  in  Hurri¬ 
cane  Allen  on  6  August  have  not  been  well  documented 
in  previous  studies,  except  that  the  vertical  velocity  data 
had  been  used  by  Jorgensen  et  al.  (1985).  Hurricane 
Allen  originated  from  a  tropical  wave  that  previously 
moved  off  the  African  coastline  on  30  July.  It  formed  as 
a  tropical  depression  on  1  August,  moving  toward  the 
west  at  nearly  20  mph.  The  storm  became  a  category  5 
storm  south  of  Puerto  Rico  on  5  August  and  stayed  at 
that  strength  for  longer  than  a  day.  It  weakened  to  a 
category  4  hurricane  on  6  August  but  regained  category 
5  intensity  again  over  a  day.  As  Allen  entered  the  Gulf 
it  weakened  back  to  a  category  4  but  once  again  quickly 
attained  category  5  status.  During  six  of  the  seven  days 
(4-10  August  1980)  when  Hurricane  Allen  was  in  the 
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Fig.  1.  Plots  of  (a)  aircraft  altitude,  (b)  temperature,  (c)  surface  wind  speed,  (d)  radial  wind  velocity,  (e)  tangential 
wind  velocity,  and  (f)  vertical  velocity  during  Hurricane  Hugo  at  1720-1729  UTC 15  Sep  1989.  The  gray  lines  denote 
the  three  time  intervals  investigated. 


Caribbean  and  Gulf  of  Mexico,  a  total  of  seven  research 
missions  were  conducted  by  the  two  NO  A  A  P3  aircraft 
(N42RF  and  N43RF).  An  overview  of  the  synoptic  con¬ 
ditions  and  a  detailed  account  of  the  field  programs  can 
be  referred  to  Jorgensen  (1984)  and  Marks  (1985). 

Figure  2  shows  the  horizontal  aircraft  track  that  is 
superimposed  on  a  composite  radar  reflectivity  image 
that  is  from  Marks  (1985).  Also  shown  in  Fig.  2  is  the 
aircraft  altitude  as  a  function  of  time  in  UTC.  During  the 
period  between  1500  to  1640  UTC,  N42RF  was  flown 
at  a  nearly  constant  altitude  of  around  450  m.  Until 
1650  UTC,  the  aircraft  remained  below  1  km  above  the 
sea  surface.  Figure  3  shows  the  variables  of  interest 
during  the  low-level  mission.  The  gray  lines  at  the  bot¬ 
tom  of  each  panel  represent  the  time  intervals  selected 
to  determine  the  scales  of  turbulent  eddies  and  turbulence 


parameters.  There  are  a  total  of  12  time  intervals  (or  flux 
runs)  selected  for  analysis,  four  of  which  are  in  the  eyewall 
region.  The  remaining  time  intervals  were  taken  out¬ 
side  the  eyewall,  and  most  of  them  correspond  to  times 
that  the  aircraft  was  between  the  outer  rainbands.  These 
remaining  time  intervals  are  used  to  compare  with  the 
CBLAST  data  at  similar  altitudes.  Note  that  all  the  time 
intervals  are  chosen  according  to  the  ogive  criterion 
discussed  in  section  3. 

From  the  radial  wind  velocity,  we  can  see  that  N42RF 
was  flown  mostly  within  the  inflow  layer  for  the  periods 
of  interest  for  both  the  Allen  and  Hugo  flights.  However, 
we  are  uncertain  if  the  flights  were  within  the  hurricane 
boundary  layer  because  there  is  no  vertical  sounding 
observation  in  both  cases.  According  to  previous  studies 
examining  vertical  profiles  of  kinematic  variables  in 
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Fig.  2.  Plots  of  (a)  aircraft  altitude  and  (b)  horizontal  aircraft  track  that  is  superimposed  on  top  of  a  composite  of 
a  radar  reflectivity  image  from  Marks  (1985),  on  6  Aug  1980  in  Hurricane  Allen. 


intense  hurricane  eyewalls  (e.g.,  Franklin  et  al.  2003; 
Powell  et  al.  2003),  it  is  likely  that  N42RF  was  near  the 
height  of  the  maximum  wind  speed.  We  assume  that 
N42RF  is  close  to  the  top  of  the  boundary  layer  for  the 
remainder  of  the  paper  (e.g.,  Kepert  and  Wang  2001; 
Bell  and  Montgomery  2008;  Smith  et  al.  2009;  Zhang 
et  al.  2009;  Smith  and  Montgomery  2010). 

Tables  1  and  2  ,  respectively,  summarize  the  measure¬ 
ments  and  calculations  for  the  time  intervals  of  the 
flights  into  Hurricanes  Hugo  and  Allen.  Overall,  the 
time- averaged  mean  wind  speeds  obtained  at  flight  level 
vary  from  7  to  64  m  s-1.  For  the  Hugo  flight,  the  av¬ 
eraged  near-surface  wind  speeds  (1/sfmr),  determined 
from  measurements  by  the  SFMR,  vary  from  21  to 
40  m  s”1.  The  stability  parameter,  z/L,  is  also  shown  in 
Tables  1  and  2,  where  z  is  the  aircraft  altitude,  and  L  is 
the  Obukhov  length  calculated  using  an  iterative  method 
(Drennan  et  al.  2007);  the  results  indicate  that  the  ther¬ 
modynamic  conditions  were  near  neutral  but  slightly 
unstable  in  the  boundary  layer.  When  estimating  L,  we 
use  the  surface  exchange  coefficients  of  momentum  and 
heat  from  Powell  et  al.  (2003)  and  Zhang  et  al.  (2008)  to 
estimate  surface  fluxes.  It  is  assumed  also  that  the  sea 
surface  temperature  is  28°C.  For  the  data  in  Hurricane 
Allen,  there  are  no  surface  wind  data,  and  we  only  have 
data  at  one  level,  so  L  was  not  estimated. 

Turbulent  fluxes  of  momentum  at  flight  level  are  cal¬ 
culated  using  the  eddy-correlation  method  for  each  time 
interval  as  follows: 


T  =  p{-w'v[  i  -  w'v'J),  (1) 

where  prime  indicates  turbulent  fluctuations;  w,  vh  and 
vr  represent  vertical,  tangential,  and  radial  component 
velocities,  respectively;  p  is  the  air  density;  and  an  over¬ 
bar  represents  a  time-averaged  operator.  TKE  is  com¬ 
puted  using  the  turbulent  fluctuations  in  the  form  of 

e=f(<  +  ^  +  ^2).  (2) 

Note  that  enthalpy  fluxes  cannot  be  estimated  in  this 
study  because  of  the  sensor  limitation  in  measuring  hu¬ 
midity  and  potential  wetting  of  the  temperature  sensor. 

Turbulent  fluctuations  are  determined  by  detrending 
the  time  series  of  the  three  wind  components  using  a 
least  squares  fitting  method.  A  high-pass  filter  with  a 
cutoff  at  0.01  Hz  was  applied  before  the  detrending. 
When  we  calculate  p,  we  use  the  temperature  measured 
by  the  Rosemont  temperature  sensor.  It  has  been  re¬ 
ported  by  Eastin  et  al.  (2002a, b)  that  there  is  usually 
a  wetting  error  in  the  temperature  data  during  eyewall 
penetrations.  We  corrected  the  wetting  error  following 
the  Eastin  et  al.  method.  We  found  that  the  influence  of 
the  wetting  error  on  the  density  calculation  is  very  small 
(—1%),  which  is  nearly  negligible. 

Because  the  vertical  eddy  diffusivity  is  a  key  param¬ 
eter  in  a  number  of  planetary  boundary  layer  (PBL)  pa¬ 
rameterization  schemes  used  in  numerical  as  well  as 
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Fig.  3.  Plots  of  (a)  aircraft  altitude,  (b)  vertical  velocity, 
during  Hurricane  Allen  at  1501-1650  UTC  6  Aug  1980 


(c)  tangential  wind  velocity,  and  (d)  radial  wind  velocity 
.  The  gray  lines  denote  the  time  intervals  investigated. 


most  theoretical  models  of  hurricanes  (e.g.,  Kepert  2001; 
Smith  et  al.  2008;  Bryan  and  Rotunno  2009b;  Smith  and 
Thomsen  2010),  we  believe  it  worthwhile  to  estimate 
the  vertical  diffusivity  here  using  the  data  in  Hurricanes 
Hugo  and  Allen.  Typically,  the  vertical  eddy  diffusivity  is 
defined  as 


K=  1^1 


-l 


(3) 


vertical  wind  shear.  For  the  case  of  Hurricane  Hugo 
when  the  SFMR  wind  data  are  available,  the  wind  shear 
is  estimated  from  the  difference  between  the  flight-level 
winds  and  the  surface  winds.  However,  for  the  Allen 
case  we  could  not  estimate  K  using  the  direct  method. 

Previous  studies  (e.g.,  Hanna  1968)  provided  an  al¬ 
ternative  method  for  estimating  K  in  the  form  of 

Kx=daw,  (4) 


where  V  is  the  mean  wind  speed,  and  z  is  the  altitude. 
We  refer  to  Eq.  (3)  as  the  direct  method  hereafter.  In 
this  method  for  estimating  K ,  one  needs  to  calculate  the 


where  aw  is  the  standard  deviation  of  the  vertical  wind 
velocity,  /  is  the  vertical  mixing  length  scale  defined 
as  l  =  (T^/s,  and  c  is  a  constant.  Here,  we  use  c  =  0.41 


Table  1.  Summary  of  data  and  calculations  for  the  three  time  intervals  (1, 2,  and  3)  in  Hurricane  Hugo  on  15  Sep  1989.  The  variables  are 
as  follows:  the  number  of  the  time  intervals  (#);  start  time  of  the  time  interval  (Ts  in  UTC  hours  and  minutes);  end  time  of  the  time  interval 
(Tnd);  mean  altitude  (z  in  m);  mean  flight-level  wind  speed  (Uz  in  m  s_1);  mean  10-m  wind  speed  (ESfmr);  stability  (where  L  is  the 
Obukhov  length);  momentum  flux  (|r|  in  N  m-2);  dissipation  rate  (e  in  cm2  s-3);  TKE  ( e  in  m2  s~2);  standard  deviation  of  the  vertical 
velocity  (aw  in  m  s_1);  vertical  mixing  length  scale  (/  in  m);  peak  horizontal  length  scale  (/phl  in  m);  and  vertical  eddy  diffusivity  (K  in 
m2  s-1).  Here,  K,  and  K2  are  the  eddy  diffusivity  estimated  using  the  direct  method,  Hanna’s  method,  and  the  TKE  closure  method, 
respectively.  Time  interval  2  is  in  the  eyewall.  Note  that  time  interval  3  is  not  a  good  run  for  flux  calculation.  The  results  for  time  interval  3 
are  given  only  for  comparison  purposes,  which  should  not  be  used  to  evaluate  model  simulation  outputs. 


# 

T 

1  s 

Tnd 

z 

Uz 

Esfmr 

z/L 

M 

e 

(jw 

8 

/ 

^phl 

K\ 

k2 

K 

1 

1720 

1724 

463 

28.2 

21.5 

-0.44 

0.34 

4.46 

0.84 

81.1 

74.5 

2490.7 

25.6 

26.3 

23.1 

2 

1725 

1728 

436 

59.0 

40.0 

-0.06 

4.70 

25.8 

2.68 

1778 

108.6 

2189.2 

118.6 

116.7 

109.4 

3 

1724 

1729 

422 

52.2 

39.1 

-0.08 

26.2 

258.2 

4.80 

5001 

221.0 

2142.8 

435.7 

497.6 

824.0 
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Table  2.  Summary  of  data  and  calculations  for  the  12  time  intervals  in  Hurricane  Allen  on  6  Aug  1980.  The  variables  are  as  follows:  #; 
Ts  (in  UTC  hours  and  minutes);  Tnd;  z  (in  m);  Uz  (in  m  s-1);  |r|  (in  N  m-2);  s  (in  cm2  s~3);  TKE  (e  in  m2  s~2);  crw  (in  m  s-1);  /  (in  m);  /phl 
(in  m);  and  K  (in  m2  s_1).  Kx  and  K2  are  the  eddy  diffusivity  estimated  using  Hanna’s  method  and  the  TKE  closure  method,  respectively. 


Note  that  time  intervals  9- 

-12  are  the  ones  in  the  eyewall. 

# 

Ts 

Tnd 

z 

Uz 

]H 

e 

(jw 

£ 

/ 

^phl 

K\ 

k2 

1 

1505 

1508 

482.6 

24.5 

0.28 

1.65 

0.78 

76.0 

63.6 

900.3 

24.9 

10.8 

2 

1508 

1512 

442.6 

25.1 

0.53 

3.36 

0.96 

132.8 

66.0 

911.4 

47.8 

25.5 

3 

1536 

1539 

473.0 

18.6 

0.13 

1.75 

0.70 

34.4 

98.1 

888.2 

11.5 

26.7 

4 

1542 

1545 

497.8 

13.3 

0.11 

1.32 

0.48 

26.9 

40.7 

1812.1 

10.2 

19.4 

5 

1546 

1549 

498.5 

11.0 

0.21 

1.39 

0.54 

56.3 

28.7 

1904.7 

18.6 

10.3 

6 

1550 

1553 

498.3 

7.2 

0.12 

1.57 

0.63 

132.3 

19.0 

3040.5 

10.8 

5.6 

7 

1613 

1616 

422.1 

19.3 

0.41 

3.37 

0.93 

64.4 

125.2 

859.9 

18.3 

52.8 

8 

1648 

1652 

847.7 

40.1 

0.88 

4.06 

1.18 

94.6 

86.9 

2947.6 

20.5 

26.1 

9 

1517 

1520 

484.6 

39.7 

2.74 

13.1 

2.20 

1243 

85.7 

951.3 

77.3 

41.7 

10 

1528 

1533 

631.1 

37.9 

0.99 

10.7 

1.66 

399 

114.6 

2521.2 

78.0 

86.7 

11 

1623 

1627 

448.7 

40.9 

1.76 

15.7 

2.46 

1191 

125.0 

1760.1 

126.1 

62.1 

12 

1632 

1637 

844.9 

64.2 

2.89 

20.0 

2.38 

1900 

71.0 

763.1 

69.2 

62.4 

following  Nieuwstadt  (1984).  Equation  (4)  will  be 
hereafter  referred  to  as  Hanna’s  method.1  Although 
Hanna’s  method  is  well  suited  for  boundary  layers 
where  similarity  theory  is  a  plausible  approximation, 
this  method  has  been  used  successfully  in  other  condi¬ 
tions,  such  as  in  the  boundary  layer  over  the  forests 
where  conditions  are  not  homogeneous  (Lee  1996).  In 
this  study,  we  will  assume  that  Hanna’s  method  can  be 
used  also  in  the  low-level  troposphere  in  hurricane 
conditions.  The  limitations  of  our  estimates  are  dis¬ 
cussed  in  sections  5  and  6. 

The  rate  of  dissipation  s  is  usually  estimated  from  the 
velocity  spectra  in  the  form  of 

e  =  a-3/2^[/5(/)]3/2;  (5) 

where  /  is  the  frequency,  a  is  the  one-dimensional 
Kolmogorov  constant,  U  is  the  true  airspeed  relative  to  the 
aircraft,  and  S  is  the  power  spectral  density  of  the  along- 
wind  component  wind  velocity.  In  this  study,  we  conser¬ 
vatively  use  a  =  0.5  following  the  laboratory  guidance 
provided  by  Sreenivasan  (1995).  In  the  estimation  of  e, 
we  assume  the  presence  of  an  inertial  subrange  in  the 
power  spectrum. 

In  the  PBL  parameterization  schemes  used  in  numer¬ 
ical  models,  especially  in  the  TKE  closure-type  schemes 
(Holt  and  Raman  1988),  eddy  diffusivity  is  usually  re¬ 
lated  to  TKE  and  the  rate  of  dissipation,  in  the  form  of 

K2  =  C2  e2/e,  (6) 


1  This  is  expected  also  from  basic  turbulence  phenomenology 

(e.g.,  Frisch  1996). 


where  c2  is  an  empirical  parameter.  Here,  we  use  c2  = 
0.03,  following  Detering  and  Etling  (1985),  Beljaars  et  al. 
(1987),  and  Lee  (1996).  Again,  in  Eq.  (6)  the  estimation 
of  K  involves  the  estimation  of  the  rate  of  dissipation, 
which  we  believe  to  have  the  largest  uncertainty  among 
all  of  the  above-mentioned  parameters.  The  uncertainty 
involved  in  the  estimation  of  these  turbulence  parameters 
is  discussed  in  sections  3  and  5. 

3.  Quality  control  and  potential  errors 
in  the  data  analysis 

The  most  important  issue  to  note  is  that  the  flight- 
level  data  used  in  this  study  have  a  sampling  rate  of  1  Hz. 
Although  these  data  were  used  widely  in  previous  stud¬ 
ies  to  understand  various  aspects  of  hurricane  dynamics, 
thermodynamics,  and  vortex  structure  (e.g.,  Willoughby 
et  al.  1982;  Jorgensen  1984;  Willoughby  1990;  Kossin  and 
Eastin  2001;  Eastin  et  al.  2005a, b;  Mallen  et  al.  2005),  they 
have  been  rarely  used  to  estimate  turbulence  properties. 
One  potential  major  concern  with  the  1-Hz  data  for 
turbulence  studies  is  the  sampling  problem — that  is,  the 
1-Hz  data  cannot  ordinarily  capture  the  whole  spectrum 
of  the  turbulent  energy.  This  limitation  notwithstanding, 
we  believe  that  the  1-Hz  data  such  as  those  collected  in 
Hurricanes  Hugo  and  Allen  still  have  their  merit  for  es¬ 
timating  turbulence  parameters.  Our  belief  is  supported 
by  a  detailed  comparison  of  the  momentum  fluxes,  vari¬ 
ances,  and  TKEs  estimated  using  the  1-  and  40-Hz  data 
for  identical  flux  runs  during  the  CBLAST  experiment. 

As  an  example,  Fig.  4  shows  a  comparison  between 
the  wind  data  with  1-  and  40-Hz  sampling  frequencies. 
The  data  shown  in  Fig.  4  are  from  a  flux  run  that  was 
collected  by  NO  A  A  N43RF  aircraft  in  Hurricane  Frances 
at  1920  UTC  1  September  2004.  The  measurements  were 
taken  between  the  outer  rainbands  of  Hurricane  Frances 
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Fig.  4.  Comparison  of  the  40-  (black)  and  1-Hz  (blue)  wind  data,  (a)  Time  series  comparison  from  a  typical  flux  run 
at  460  m  on  1  Sep  2004  in  Hurricane  Frances,  (b)  Frequency  spectra  from  the  flux  run  of  (a),  (c)  Frequency  cospectra 
of  uw.  (d)  Cumulative  sum  or  ogive  of  the  spectra,  (e)  Cumulative  sum  or  ogive  of  cospectra. 


as  part  of  the  CBLAST  mission.  The  mean  flight-level 
wind  speed  of  this  flux  run  is  32  m  s_1.  The  flux  data  for 
this  flight  has  been  presented  elsewhere  by  Zhang  et  al. 
(2009).  We  note  that  this  flight  is  the  only  one  that  has 
measurements  between  400  and  500  m  during  CBLAST. 

The  time  series  of  the  along-wind  component  wind 
velocity  u  in  1  and  40  Hz  are  compared  in  Fig.  4a;  the 
analyzed  data  show  good  agreement.  The  spectra  and 
cospectra  comparisons  are  shown  next  in  Figs.  4b  and  4c, 
respectively.  We  see  that  there  is  a  good  agreement  in 
the  spectra  and  cospectra  between  the  1-  and  40-Hz  data 
at  frequencies  <0.5  Hz,  the  equivalent  of  the  Nyquist 
frequency  for  1-Hz  data.  It  is  noticed  that  there  is  a  ten¬ 
dency  for  the  1-Hz  data  to  decay  faster  at  a  frequency 
range  between  0.4  and  0.5  Hz.  Beyond  the  Nyquist  fre¬ 
quency,  however,  the  1-Hz  data  is  unable  to  measure 
turbulent  eddies  smaller  than  around  200  m  in  space. 
Now,  to  estimate  the  spatial  scales  of  the  turbulent 
eddies,  we  first  applied  Taylor’s  frozen  hypothesis 
to  transfer  the  frequency  to  the  wavenumber  domain, 


defined  as  the  frequency  divided  by  the  true  airspeed 
relative  to  the  aircraft,  which  is  on  the  order  of  100- 
150  m  s-1.  The  scale  is  then  the  reciprocal  of  the  wave- 
number. 

The  cumulative  spectral  and  cospectral  sums,  or 
“ogives,”  are  shown  in  Figs.  4d  and  4e,  respectively.  The 
ogive  plots  have  been  widely  used  in  turbulence  studies 
to  determine  useful  time  intervals  suitable  for  turbu¬ 
lent  flux  calculations  (e.g.,  Zhang  et  al.  2009).  The  use 
of  ogive  curves  when  estimating  fluxes  provides  insight 
into  turbulent  structure  and  spatial  scales  that  contrib¬ 
ute  to  the  turbulent  transport.  The  flatness  of  the  ogive 
curve  at  low-  and  high-end  frequencies  indicates  that 
the  energy -containing  scale  of  the  variance  and  flux  pro¬ 
cesses  are  well  sampled.  If  the  ogive  curve  approaches 
asymptotically  a  single  value,  then  homogeneity  or  sta- 
tionarity  of  the  flow  is  met,  and  the  asymptotic  value 
of  the  ogive  represents  the  total  variance  and  covari¬ 
ance  or  flux  of  momentum.  Upon  examining  Figs.  4d 
and  4e,  we  see  that  the  total  variance  and  covariance  are 
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Table  3.  Summary  of  data  and  calculations  for  the  10  time  intervals  in  Hurricane  Frances  on  1  Sep  2004.  The  variables  are  as  follows:  #; 
Ts  (in  UTC  hours  and  minutes);  Tnd;  z  (in  m);  Uz  (in  m  s-1);  |r|  (in  N  m-2);  e  (in  cm2  s“3);  TKE  ( e  in  m2  s~2);  crw  (in  m  s-1);  /  (in  m);  and  K 
(in  m2  s_1).  The  subscript  1  represents  1-Hz  data,  and  subscript  40  represents  40-Hz  data.  K±  and  K2  are  the  eddy  diffusivity  estimated 
using  Hanna’s  method  and  the  TKE  closure  method,  respectively.  Vertical  mixing  length  and  Ks  are  calculated  using  the  40-Hz  data. 


# 

Ts 

Tnd 

z 

Uzl 

UZ40 

e\ 

^40 

Mi 

M« 

Owl 

O\v40 

Sl 

S40 

/ 

K\ 

k2 

1 

1743 

1746 

484.1 

20.6 

20.3 

0.83 

1.14 

0.05 

0.08 

0.21 

0.48 

7.8 

25.1 

44.7 

7.6 

15.5 

2 

1747 

1750 

477.6 

20.6 

20.3 

0.77 

1.25 

0.15 

0.19 

0.23 

0.65 

11.2 

25.5 

106.3 

16.5 

18.4 

3 

1751 

1754 

484.5 

21.0 

20.4 

0.79 

1.36 

0.21 

0.24 

0.38 

0.80 

11.6 

30.6 

144.2 

42.2 

15.6 

4 

1755 

1759 

484.1 

22.2 

21.6 

1.19 

1.81 

0.27 

0.44 

0.40 

0.81 

11.6 

15.6 

336.3 

39.4 

63.0 

5 

1801 

1805 

451.2 

21.2 

21.7 

1.16 

1.78 

0.19 

0.23 

0.31 

0.87 

39.5 

65.4 

102.1 

20.4 

14.5 

6 

1907 

1910 

563.2 

29.4 

28.8 

2.41 

3.33 

0.53 

0.65 

0.61 

1.03 

44.2 

97.6 

113.2 

57.8 

34.1 

7 

1911 

1913 

529.7 

33.4 

32.7 

1.28 

1.71 

0.18 

0.26 

0.33 

0.67 

8.0 

11.7 

253.5 

22.8 

74.5 

8 

1914 

1917 

498.4 

33.2 

32.6 

1.04 

1.57 

0.22 

0.24 

0.36 

0.80 

17.0 

29.1 

173.4 

21.3 

25.4 

9 

1917 

1921 

486.7 

32.9 

32.2 

1.86 

2.43 

0.39 

0.51 

0.56 

0.92 

28.7 

52.3 

147.5 

45.7 

33.8 

10 

1925 

1928 

467.0 

33.3 

33.7 

2.18 

4.05 

0.31 

0.47 

0.42 

1.07 

133.0 

231.9 

52.1 

42.4 

21.3 

comparable  to  each  other  for  the  1-  and  40-Hz  data.  It  is 
evident  also  that  the  1-Hz  data  somewhat  underestimates 
the  variance  and  covariance. 

On  average,  we  found  that  the  1-Hz  Frances  data 
capture  approximately  75%  of  the  total  momentum  flux 
and  70%  of  the  TKE  (Table  3).  In  the  analysis  of  the 
Allen  and  Hugo  data,  an  empirical  correction  is  applied 
to  the  flux  and  TKE  calculations  to  compensate  for  the 
bias  determined  from  the  Frances  data.  We  have  assumed 
that  the  turbulence  characteristics  at  the  similar  vertical 
levels  and  locations  in  Frances,  Hugo,  and  Allen  behave 
in  a  similar  manner.  We  recognize  that  there  is  uncer¬ 
tainty  in  the  correction,  especially  in  the  eyewall  region 
where  the  40-Hz  data  are  unavailable,  but  this  approach 
provides  our  best  estimates. 

Boundary  layer  turbulence  studies  usually  assume  that 
turbulent  properties  such  as  variances  and  fluxes  are  en¬ 
semble  averages  such  that  the  average  of  a  number  of 
individual  samples  of  data  is  under  identical  conditions. 
In  aircraft  measurements,  it  is  generally  assumed  also 
that  time  or  space  averages  converge  to  ensemble  av¬ 
erages  when  the  averaging  time  becomes  sufficiently 
long  or  the  averaging  length  becomes  sufficiently  large 
(e.g.,  Wyngaard  1986;  Lenschow  1986).  In  the  present 
case,  the  time  intervals  suitable  for  analysis  are  relatively 
short;  these  intervals  were  selected  based  on  the  quality 
control  through  the  ogives  criterion  discussed  above  (fur¬ 
ther  details  are  given  in  section  4).  The  short  time  in¬ 
terval  used  in  the  calculations  can  yield  uncertainty  of 
variance  and  covariance  fluxes  according  to  Mann  and 
Lenschow  (1994)  and  Mahrt  (1998).  We  present  the  de¬ 
tailed  error  analysis  in  section  5  when  we  present  the 
results  of  the  estimated  turbulence  parameters. 

4.  Spectral  analysis  of  eddy  structure 

The  spectra  of  the  three-dimensional  wind  velocities 
for  the  three  time  intervals  in  Hurricane  Hugo  are  plotted 


in  Fig.  5.  Note  that  Marks  et  al.  (2008)  conducted  spec¬ 
tral  analysis  for  the  whole  leg  that  contains  the  three 
time  intervals.  Here,  a  more  detailed  spectral  analysis 
for  each  of  the  three  time  intervals  is  shown  separately. 
Figure  5  indicates  that  the  spectra  of  the  three-component 
wind  velocity  for  all  the  three  time  intervals  tends  to 
follow  the  —  5/3  Kolmogorov  law  at  frequencies  >0.2  Hz. 
Because  the  time  intervals  of  the  eyewall  penetrations 
are  in  the  crosswind  direction,  we  cannot  determine 
if  there  is  a  ratio  of  %  between  the  cross-stream  and 
along-stream  velocity  spectra  (Zhang  2010).  To  estimate 
the  turbulent  dissipation  rate  using  Eq.  (4),  we  must 
assume  that  the  inertial  subrange  exists.  It  is  evident  in 
Fig.  5  that  the  peak  energy  in  the  wind  velocity  spectra 
increases  somewhat  from  time  intervals  1  to  3.  This  in¬ 
crease  of  energy  in  the  spectra  is  consistent  with  the 
increase  of  the  rate  of  dissipation  from  intervals  1  to  3, 
as  shown  in  Table  1. 

The  cospectra  of  the  vt  and  vr  components  of  the  mo¬ 
mentum  fluxes  for  time  intervals  1  and  2  of  Hurricane 
Hugo  flight  are  shown  next  in  Figs.  6a  and  6b,  respec¬ 
tively.  The  corresponding  ogives  of  the  cospectra  are 
shown  in  Figs.  6c  and  6d.  To  estimate  the  spatial  scale 
of  turbulent  eddies,  we  plot  the  cospectra  and  ogives 
in  the  wavenumber  domain.  Since  the  ogive  curves  are 
nearly  flat  at  low  and  high  wavenumbers,  this  supports 
the  assumption  of  the  stationarity  for  time  intervals  1 
and  2.  The  dominant  peaks  in  the  cospectral  plots  are 
associated  with  turbulent  eddies  that  contain  most  of  the 
momentum  fluxes.  As  mentioned  earlier,  the  scales  of 
the  dominant  eddies  can  be  estimated  as  the  reciprocal  of 
the  wavenumbers.  During  time  interval  1  (Fig.  6a),  for 
example,  the  peaks  in  the  vtw  cospectral  plots  at  wave- 
numbers  4  X  10-4  m-1  and  6  X  ICE4  m-1  correspond 
to  length  scales  of  approximately  2.5  and  1.4  km,  re¬ 
spectively.  During  time  interval  2  (Fig.  6b,  the  eyewall 
penetration  leg),  the  dominant  peaks  of  the  vtw  cospectra 
at  wavenumbers  4.5  X  ICE4  m-1  and  1.5  X  ICE3  m-1 


May  2011 


ZHANG  ET  AL. 


1455 


time  interval  1  time  interval  2 


time  interval  3 


Fig.  5.  Spectra  of  three  velocity  components  for  time  intervals  (a)  1,  (b)  2,  and  (c)  3  in  Hurricane  Hugo.  The  dotted 

line  shows  a  —  5/3  slope. 


correspond  to  length  scales  of  approximately  2.2  km  and 
670  m,  respectively. 

The  cospectra  and  ogives  plots  for  the  four  eyewall 
penetrations  legs  in  Hurricane  Allen  are  shown  in  Fig.  7. 
It  is  evident  from  the  plotted  data  that  the  eddy  struc¬ 
tures  in  those  legs  are  comparable  to  that  found  in  the 
eyewall  penetration  leg  of  Hurricane  Hugo.  The  dom¬ 
inant  peaks  of  the  vtw  and  vrw  cospectra  take  place 
at  scales  of  approximately  500-3000  m.  From  the  cor¬ 
responding  ogive  curves  it  appears  also  that  eddies  with 
such  scales  contribute  most  to  the  total  momentum  fluxes. 
Tables  1  and  2  list  the  peak  horizontal  length  scales  of 
dominant  eddy  for  all  the  time  intervals.  The  results 
suggest  that  there  is  no  apparent  dependence  of  the  peak 
horizontal  length  scale  on  wind  speed. 

Figure  8  shows  the  cospectra  of  the  momentum  fluxes 
and  the  corresponding  ogives  for  time  interval  3  of  the 
Hugo  flight.  From  the  ogive  curves  in  Fig.  8b,  the  con¬ 
tribution  of  turbulent  eddies  to  the  total  flux  in  the  wave- 
number  range  1.9  X  10-4  m-1  to  3  X  10-3  m_1  increases 
as  the  wavenumber  decreases.  However,  there  is  a  sharp 


increase  in  the  contribution  to  the  total  energy  for  wave- 
number  <1.9  X  10_4m_1,  indicating  that  the  vertical 
momentum  transport  is  strongly  influenced  by  the  large- 
scale  features,  possibly  the  convective-scale  eddies  with 
scales  >5  km.  The  inclusion  of  the  EVM  increases  the 
total  momentum  flux  by  almost  an  order  of  magnitude 
compared  to  that  of  time  interval  2.  On  the  other  hand, 
because  of  the  sharp  increase  in  the  ogive  plots,  Fig.  8 
indicates  that  time  interval  3  is  not  a  good  run  for  the  flux 
calculation.  Based  upon  the  discussion  and  interpreta¬ 
tion  presented  in  Marks  et  al.  (2008),  the  EVMs  appear 
to  behave  in  some  respects  like  quasi-two-dimensional 
vortex  structures  that  are  embedded  within  the  three- 
dimensional  turbulence  on  the  inside  edge  of  the  eyewall. 
Such  features  are  believed  to  have  a  transient  nature  but 
possess  very  intense  signatures  locally  (Schubert  et  al. 
1999;  Kossin  and  Schubert  2001;  Montgomery  et  al. 
2002;  Yau  et  al.  2004;  Montgomery  et  al.  2006).  If  the 
EVMs  were  included  in  the  turbulence  analysis,  the  ho¬ 
mogeneity  or  stationarity  condition  would  not  be  satis¬ 
fied,  indicating  that  it  is  unwise  to  include  EVM  in  the 
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Time  interval  1 


Time  interval  2 


wavenumber  [rrf1] 

Fig.  6.  Cospectra  of  the  vertical  velocity  w  with  the  tangential  vt  and  radial  vr  velocity  and  the  corresponding 
cumulative  sums  or  ogives  for  time  intervals  1  and  2,  respectively,  in  Hurricane  Hugo. 


flux  calculation.  We  point  out  here  that  we  have  listed 
the  value  of  the  vertical  momentum  flux  for  time  interval 
3  in  Table  1  for  comparison  purposes  only.  We  think  that 
time  interval  2  is  more  like  a  typical  eyewall  boundary 
layer  leg  at  450  m  than  time  interval  3.  Thus,  the  numerical- 
modeling  community  should  be  cautious  when  using  the 
results  of  time  interval  3  of  the  Hugo  flight  to  compare 
with  model  output. 

5.  Estimation  of  turbulence  parameters 
and  error  analysis 

The  calculations  of  TKE  and  momentum  flux  using 
Eqs.  (1)  and  (2)  are  summarized  in  Tables  1  and  2,  re¬ 
spectively,  for  the  Hugo  and  Allen  data.  The  data  outside 
the  eye  walls  of  Hurricanes  Hugo  and  Allen  show  good 
agreement  with  the  40-Hz  Frances  data.  This  agreement 
suggests  that  the  method  employed  for  correcting  the 
1-Hz  data  is  sound.  The  momentum  flux  for  the  eye- 
wall  penetration  leg  in  Hurricane  Hugo  is  approxi¬ 
mately  4.7  N  m-2  with  a  flight-level  mean  wind  speed 
of  59  m  s_1.  The  mean  momentum  flux  of  the  two  eyewall 


legs  at  around  450  m  for  Hurricane  Allen,  with  a  flight- 
level  mean  wind  speed  of  approximately  40  m  s_1,  is 
2.3  N  m-2.  Overall,  the  momentum  fluxes  in  the  eyewall 
legs  are  generally  5-10  times  those  found  in  the  outer- 
core  runs.  The  TKE  values  for  the  eyewall  legs  fall  be¬ 
tween  10  and  25  m2  s-2;  these  values  are  consistent  with 
the  independently  derived  results  of  Lorsolo  et  al.  (2010) 
who  used  Doppler  radar  data  to  map  TKE  in  several 
hurricanes.  The  TKE  in  the  eyewall  region  obtained 
using  the  flight-level  data  is  roughly  7  times  that  in  the 
legs  outside  the  eyewall;  this  is  also  consistent  with  the 
Lorsolo  et  al.  study. 

Figures  9a  and  9b  show  the  TKE  and  momentum 
fluxes  as  a  function  of  the  flight-level  mean  wind  speed, 
respectively.  Also  shown  are  the  values  of  TKE  and  mo¬ 
mentum  fluxes  determined  from  the  40-Hz  data  ob¬ 
tained  in  Hurricane  Frances  (2004).  Based  upon  these 
results,  it  appears  that  the  momentum  flux  and  TKE 
increase  with  the  increasing  flight-level  wind  speed,  con¬ 
firming  that  the  turbulent  flux  and  energy  are  strongest 
in  the  eyewall  region. 

As  mentioned  in  section  3,  there  are  sources  of  error 
that  are  involved  in  the  flux  estimation.  Following  the 


May  2011 


ZHANG  ET  AL. 


1457 


10  10  10 
wave  number  [m'1] 

Fig.  7.  As  in  Fig.  6,  but  for  the  (left  to  right)  four  time  intervals  9, 10, 11,  and  12  during  the  eyewall  penetrations  of  Hurricane  Allen. 


method  described  by  Drennan  et  al.  (2007),  the  sampling 
error  is  given  by 

aF/F  =  aFzll2U-ll2T-112,  (7) 

where  aF  is  the  standard  deviation  of  flux  estimates,  aF 
is  a  constant,  Ua  is  the  true  airspeed  relative  to  the  air¬ 
plane,  T  is  the  sampling  interval  (s),  and  z  is  the  altitude 
(m).  For  momentum  flux,  Sreenivansan  et  al.  (1978) 
estimated  aF  =  3.  The  measurements  in  this  study  were 
carried  out  at  altitudes  between  430  and  800  m,  with 
a  true  airspeed  of  Ua  ~  120  ms-1  and  a  short  duration 
of  T-  180  s.  Using  Eq.  (7),  the  expected  variability  of 
the  flux  estimates  using  the  Hugo  and  Allen  data  is  then 
85%.  This  variability  is  higher  than  the  measured  vari¬ 
ability  of  53%  for  the  five  eyewall  legs  and  74%  for  the 
legs  outside  the  eyewall,  suggesting  that  the  analysis 
method  used  here  is  reasonable. 

Two  types  of  errors  arise  when  estimating  fluxes 
by  the  eddy-correlation  method:  the  systematic  error 
(erS),  which  is  linked  to  the  loss  due  to  high-pass  filter¬ 
ing;  and  random  error  (erR),  which  is  due  to  the  fact  that 
a  flight  leg  is  a  finite  sample  of  a  random  process.  We 
calculate  the  systematic  error  for  the  flux  estimation 
following  the  methodology  described  by  Mann  and 
Lenschow  (1994), 


erS  =  (F  -  Ff)/F,  (8) 

where  Ff  is  the  flux  after  the  high-pass  filter.  The  sys¬ 
tematic  error  is  found  to  be  28%,  which  is  normal  for 
aircraft  observations  especially  at  these  altitudes.  We 
estimate  the  random  error  following  Vickers  and  Mahrt 
(1997), 


Fig.  8.  As  in  Fig  6,  but  for  time  interval  3  in  Hurricane  Hugo. 
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Fig.  9.  Plots  of  TKE  and  momentum  flux  as  a  function  of  the  mean  wind  speed  at  the  flight  level  in  Hurricanes 

Allen  (X),  Hugo  (+),  and  Frances  (O). 


erR  =  ctf/F/Vn,  (9) 

where  N  is  the  number  of  observations.  The  random 
error  is  found  to  be  24%  for  both  the  eyewall  legs  and 
those  outside  the  eyewall,  which  is  in  agreement  with 
values  found  in  the  literature  (e.g.,  Mann  and  Lenschow 
1994;  Bernard-Trottolo  et  al.  2003).  Because  all  the  time 
intervals  or  flux  runs  were  thoroughly  checked  using 
ogive  criterion  and  the  cumulative  sum  analysis  method 
as  mentioned  in  section  3,  all  of  the  low-frequency  scales 
are  captured.  We  have  also  corrected  the  missing  high- 
frequency  part  of  the  energy  based  on  the  CBLAST  data. 
At  the  end,  the  overall  uncertainty  of  the  estimated  TKE 
and  momentum  flux  is  thought  to  be  around  30%. 

Using  Hanna’s  method,  we  find  that  the  vertical  mix¬ 
ing  length  scale  near  the  eyewall  region  is  approximately 
70-125  m.  The  mixing  length  estimated  using  the  Frances 
data  varies  from  50-330  m.  Similar  values  of  mixing 
length  are  observed  in  the  Allen  and  Hugo  data  in  the 
outer-core  region.  In  light  of  the  large  scatter  of  the 
data,  we  cannot  deduce  a  wind  speed  dependence  of 
the  mixing  length. 

The  vertical  eddy  diffusivity  is  estimated  using  the 
methods  mentioned  in  section  2.  Using  the  Hugo  data, 
the  eddy  diffusivity  estimated  based  on  the  three  dif¬ 
ferent  methods  are  in  reasonable  agreement  (Table  1). 


Using  the  Allen  data,  the  Ks  estimated  using  Hanna’s 
method  and  the  TKE  closure  method  also  agree  with  each 
other  (Table  2). 

In  estimating  K  using  the  direct  method  [Eq.  (3)],  we 
estimated  the  vertical  shear  using  the  wind  data  at  two 
levels.  It  is  likely  that  we  have  overestimated  the  shear, 
given  that  the  observation  altitude  may  be  close  to  the 
wind  maximum  where  the  shear  is  relatively  low  ac¬ 
cording  to  the  eyewall  dropsonde  data  shown  by 
Franklin  et  al.  (2003)  and  Powell  et  al.  (2003).  However, 
the  data  we  used  are  outside  the  radius  of  the  maximum 
wind  speed,  where  the  shear  at  450  m  is  usually  higher 
than  that  in  the  inner  eyewall  from  previous  case  studies 
(e.g.,  Kepert  2006a, b;  Schwendike  and  Kepert  2008).  We 
believe  that  our  estimation  of  the  shear  in  Hurricane 
Hugo  is  realistic. 

Figure  10  shows  the  Ks  as  a  function  of  the  mean 
flight-level  wind  speed,  using  the  Hugo  and  Allen  data  as 
well  as  the  40-Hz  Frances  data.  The  Ks  estimated  only 
using  the  Frances  data  based  on  Hanna’s  method  and 
the  TKE  closure  methods  show  good  agreement.  We 
found  that  the  Allen  data  and  Hugo  data  in  the  outer- 
core  region  are  close  to  the  Frances  data,  giving  some 
assurance  of  the  validity  of  the  bias  correction.  Consid¬ 
ering  all  the  data  investigated  in  this  work,  it  is  evident 
that  Ks  in  the  eyewall  regions  are  much  larger  than  that 
those  at  the  outer-core  regions.  For  the  eyewall  legs, 
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Fig.  10.  Plots  of  the  eddy  diffusivity  as  a  function  of  mean  wind 
speed  at  the  flight  level  for  all  the  good  runs  in  Hurricanes  Allen, 
Hugo,  and  Frances.  Here,  Ks  are  estimated  using  three  methods: 
the  direct  method  (O),  Hanna’s  method  (X),  and  the  TKE  closure 
method  (+).  The  Frances  data  are  in  blue.  The  Hugo  and  Allen 
data  are  in  the  other  color  for  different  methods. 

the  Ks  vary  between  70  and  130  m2  s_1  using  Hanna’s 
method,  and  they  vary  between  40  and  90  m2  s-1  using 
the  TKE  closure  method.  Overall,  it  is  found  that  the 
eddy  diffusivity  tends  to  increase  logarithmically  with  the 
increasing  mean  flight-level  wind  speed. 

When  using  either  Hanna’s  method  or  the  TKE  clo¬ 
sure  method,  a  large  error  arises  in  the  estimation  of 
the  rate  of  dissipation.  To  reduce  the  uncertainty  in  the  s 
estimation,  we  applied  a  correction  to  the  Hugo  and 
Allen  1-Hz  data  based  on  the  Frances  40-Hz  data.  After 
the  correction,  the  Hugo  and  Allen  data  at  the  outer- 
core  region  agree  with  the  Frances  data.  Given  that  there 
are  no  high-resolution  (>10  Hz)  data  available,  it  is  not 
possible  to  precisely  ascertain  the  uncertainty  in  s  in  the 
eyewall  region.  The  authors  are  unaware  of  any  other 
method  for  estimating  K  that  does  not  rely  on  an  esti¬ 
mation  of  e,  given  the  limited  data  used  here.  The  uncer¬ 
tainty  in  the  estimation  of  K  using  Hanna’s  method  and 
the  TKE  closure  method  also  comes  from  the  empirical 
constants.  For  instance,  we  used  c  =  0.41  following 
Nieuwstadt  (1984)  for  stable  boundary  layer,  while  the 
boundary  layer  we  studied  are  nearly  neutral.  The  overall 
uncertainty  of  the  estimated  K  using  Hanna’s  method  and 
the  TKE  closure  methods  is  thought  to  be  around  50%. 

6.  Discussion  and  conclusions 

In  this  study,  the  turbulence  characteristics  in  the  low- 
level  troposphere  (—450  m  altitude)  of  Hurricanes  Hugo 


(1989)  and  Allen  (1980)  have  been  investigated.  Tur¬ 
bulent  fluxes  of  momentum  and  TKE  were  estimated 
before  and  during  the  eyewall  penetration.  Momentum 
fluxes  and  TKE  estimated  for  the  eyewall  penetration 
leg  are  found  to  be  nearly  an  order  of  magnitude  larger 
than  those  estimated  for  the  legs  outside  the  eyewall  at 
the  same  level.  The  TKE  and  momentum  fluxes  are 
found  to  increase  with  increasing  mean  wind  speed  at 
the  same  level.  The  vertical  mixing  length  scale  is  found 
to  be  approximately  100  m  in  the  eyewall  region,  with 
slightly  smaller  values  outside  the  eyewall. 

Through  spectral  analysis,  the  spatial  scales  of  domi¬ 
nant  turbulent  eddies  in  the  eyewall  penetration  legs  are 
found  to  lie  between  500  and  3000  m.  The  turbulence  in 
the  hurricane  boundary  layer  is  three-dimensional  from 
the  spectral  analysis.  Our  analyses  indicate  that  it  is  un¬ 
wise  to  include  EVM  in  the  turbulence  parameter  esti¬ 
mation. 

We  also  estimated  the  vertical  eddy  diffusivity  K  using 
three  different  methods:  the  first  uses  the  definition  of  K, 
the  second  uses  a  theoretical  method  given  by  Hanna 
(1968),  and  the  third  uses  a  TKE  closure  method.  Us¬ 
ing  these  three  methods,  the  estimated  Ks  are  generally 
consistent  with  one  another  using  the  Hugo  data.  The 
estimated  Ks  using  the  later  two  methods  for  the  Allen 
data  also  agree  with  each  other.  Based  upon  the  rea¬ 
sonable  agreement  on  the  estimated  Ks  using  the  Hugo 
and  Allen  data  outside  the  eyewall  with  the  CBLAST 
data  in  Hurricane  Frances,  we  think  our  methodology 
for  estimating  K  using  the  available  data  is  sound.  In  the 
intense  eyewall  region,  we  found  that  K  varies  approx¬ 
imately  from  40  to  130  m2  s-1. 

Foster  (2009)  has  pointed  out  that  the  solutions  for 
a  similarity  model  of  the  hurricane  boundary  layer  are 
sensitive  to  the  specification  of  K.  This  solution  sensi¬ 
tivity  arises  not  only  from  the  magnitude  of  K  but  also 
the  variation  of  K  with  altitude  and  radius.  Our  anal¬ 
yses  suggest  that  K  increases  with  increasing  mean  wind 
speed  at  the  same  level.  The  use  of  a  constant  K  is  not 
expected  to  reproduce  realistic  turbulent  momentum 
fluxes  in  a  hurricane  boundary  layer. 

According  to  different  turbulence  parameterization 
methods  used  in  previous  theoretical  and  numerical  stud¬ 
ies  of  the  hurricane  boundary  layer,  the  maximum  K  varies 
from  38  to  101  m2  s-1  (e.g.,  Smith  1968;  Kepert  2001; 
Smith  2003;  Foster  2009).  This  range  is  comparable  to 
our  estimate  of  K  for  the  eyewall  legs.  On  the  other 
hand,  the  maximum  K  used  in  the  Medium-Range  Fore¬ 
cast  Model  (MRF)  boundary  layer  scheme  in  MM5  ex¬ 
ceeds  250  m2  s-1  in  the  eyewall  region  and  still  has 
values  larger  than  50  m2  s-1  at  radii  beyond  200  km 
(Braun  and  Tao  2000;  Smith  and  Thomsen  2010).  Based 
upon  our  analysis,  the  MRF  scheme  may  be  a  bit  too 
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diffusive,  while  the  rest  of  the  schemes  may  use  a  rela¬ 
tively  plausible  range  of  K.  However,  more  observations 
in  the  high-wind  regime  are  required  in  the  future  to 
evaluate  the  Ks  used  in  idealized  theoretical  and  op¬ 
erational  hurricane  numerical  models. 

It  is  certainly  true  that  our  method  to  estimate  K  is 
borrowed  from  the  standard  turbulence  community  for 
boundary  layers  in  low-wind  conditions  where  the  tur¬ 
bulence  characteristics  are  presumed  homogeneous  in 
the  horizontal  direction.  A  hurricane  boundary  layer  near 
the  eyewall  region  is  clearly  far  from  a  homogeneous 
regime!  In  particular,  when  the  transient  and  intense 
EVMs  exist,  the  flow  is  strongly  inhomogeneous.  These 
points  notwithstanding,  we  think  standard  turbulence 
phenomenology  should  still  provide  a  reasonable  esti¬ 
mate  of  the  turbulent  quantities  using  the  real  observa¬ 
tional  data.  While  an  eddy-resolving  numerical  simulation 
of  an  intense  hurricane  offers  an  alternative  approach 
for  estimating  the  turbulent  properties  that  is  free  of 
these  limitations  (Rotunno  et  al.  2009),  this  approach 
has  its  own  challenges.  We  think  our  analyses  provide 
a  useful  starting  point  for  the  evaluation  of  hurricane 
prediction  models.  Specifically,  our  results  offer  the 
opportunity  to  assess  boundary  layer  parameterization 
schemes  used  in  numerical  models  for  predicting  hurri¬ 
cane  intensity. 

One  clear  limitation  of  this  work  is  the  small  sample  of 
the  data  used  in  the  analysis.  In  searching  for  the  low- 
altitude  flight-level  database  collected  by  the  Hurricane 
Research  Division  during  the  last  several  decades,  the 
Allen  and  Hugo  flight-level  datasets  are  believed  to  be 
among  the  few  available  in  situ  observations  that  were 
taken  below  500  m  during  the  eyewall  penetration  of  a 
category  4  and  5  storms.  It  is  unfortunate  for  hurricane 
science  that  we  may  not  see  such  data  like  those  collected 
in  Hurricanes  Hugo  and  Allen  in  the  near  future  because 
NOAA  has  forbidden  these  types  of  flights  with  manned 
aircraft  for  obvious  safety  reasons .  To  fully  assess  the  im¬ 
pact  of  turbulent  processes  on  the  prediction  of  hurricane 
intensification  and  maximum  intensity,  and  to  further  un¬ 
derstand  the  role  of  turbulence  in  the  eyewall  dynamics 
and  thermodynamics  of  a  hurricane,  a  dedicated  field 
program  is  recommended,  possibly  with  unmanned  plat¬ 
forms  employing  advanced  turbulent  sensors  on  board. 
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